J/ip transport in QGP and p t distribution at SPS and RHIC 
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Combining the hydrodynamic equations for the QGP evolution and the transport equation for the 
primordially produced J/ip in the QGP, we investigate the J/ip transverse momentum distribution 
as well as its suppression in the i/s=17.3A GeV Pb-Pb collisions at SPS and -/s=200A GeV Au-Au 
collisions at RHIC. The two sets of equations are connected by the J/ip anomalous suppression 
induced by its inelastic scattering with gluons in the QGP. The calculated centrality dependence 
of J/ip suppression and average transverse momentum square agree well with the SPS data. From 
the comparison with the coalescence model where charm quark is fully thermalized, our calculated 
elliptical flow of the primordially produced J/ip is much smaller. This may be helpful to differentiate 
the J/ip production mechanisms in relativistic heavy ion collisions. 

PACS numbers: 25.75.-q, 12.38.Mh, 24.85. +p 



INTRODUCTION 



Many mechanisms are proposed to explain the phe- 
nomenon of charmonium suppressionjllin relativistic 
heavy ion collisions 0- In models 0, IE IE IE IE 
IE EE fllj - the charmonium is supposed to be cre- 
ated by the hard processes in the initial state, on 
the way out, it firstly collides inelastically with spec- 
tators which leads to the normal suppression|l2l Il3j]. 
The anomalous suppression is attributed to differ- 
ent mechanisms: collisions with spectators^ , instantly 
melting in Quark-Gluon Plasma (QGP)|j|, collisions 
with gluons in QGPja, |6j, or collisions with hadronic 
comovers|E IE EE EL Recently, there appeared an- 
other kind of mo dels | lH Il5| based on the assumption 
that charmonium is created at QCD hadronization ac- 
cording to statistical law, which is the extrapolation 
of the well-tested thermal model^E] f° r light hadrons. 
There is also a two component model [T^| which mixes 
these two kinds of mechanisms. With some adjustable 
parameters in each mechanism, almost all models can 
describe the J/ip anomalous suppression data at SPS 
very well. In order to differentiate these models, we 
need more information such as open charm, Xc and ip', 
and transverse momentum (p t ) distribution. In this 
Letter, we concentrate on the observables related to 
the pt distribution of J/ip, the average transverse mo- 
mentum (pfj and the anisotropic asymmetry param- 
eter (^2) of J/ip as well as the pt integrated anoma- 
lous suppression in the heavy ion collisions at SPS and 
RHIC. 

We follow the method in Refs. [E IE to treat the 
anomalous suppression. We neglect the formation time 
and assume that the charmonium is created instanta- 
neously after the binary collisions. On the way out, 
it collides firstly with spectators and then with glu- 
ons in the QGP. The former and the later are, respec- 
tively the origin of normal and anomalous suppres- 
sion in our approach. The QGP evolution is calcu- 
lated with a 2 + 1 dimensional boost invariant rela- 
tivistic hydrodynamics 19] . Since J/ip is heavy enough, 
a classical Boltzmann-type transport equation in the 
transverse phase space [l8| is used to describe the evo- 
lution of its transverse distribution function. We as- 
sume that the local equilibrium is reached at a proper 



time To when the normal suppression has ceased and 
the anomalous suppression starts. Therefore, the ini- 
tial condition of the transport equation is determined 
by the normal suppression. We neglect the elastic 
collisions between charmonium and particles in the 
medium, for the much larger charmonium mass than 
the typical temperature of the medium. 

The paper is organized as follows. We describe the 
evolution of the medium in Section 2, give the details 
of the cross sections between charmonium and particles 
in Section 3, discuss the transverse transport equation 
in Section 4, and show the numerical results in Section 
5. Finally we conclude in Section 6. 



II. MEDIUM EVOLUTION 

As in Ref.[2(|, we assume that the produced medium 
reaches local equilibrium at a proper time tq = 0.8 fm/c 
for y/s = I7.3A GeV Pb-Pb collisions at SPS and 0.6 
fm/c for y/s = 200 A GeV Au-Au collisions at RHIC. 
The consequent evolution is described with relativistic 
hydrodynamics, 



0, 



0, 



(1) 



where T^ v = (e + p)u^u v — g^ v p is the energy- 
momentum tensor and — nu^ the baryon current 
with four-velocity u M of the fluid cell, energy density 
e, pressure p and baryon density n. In our calculation, 
we use Bjorken's hydrodynamical model^El- With the 
Hubble-like longitudinal expansion and boost invari- 
ant initial condition, the hydrodynamical quantities are 
functions of the proper time r = \/t 2 — z 2 and trans- 
verse coordinates only, and the equations can be sim- 
plified as 



d T E + V • M 
d T M x + V • (M x v) 



-(E + p)/r , 
-M x /r - d x p , 
d T M y + V • (M y v) = -M v /t - d y p , 
d T R + V • (Rv) = -R/t 



(2) 



with the definitions E = (e + p)~f 2 — p, M = (e +p)j 2 v 
and R = 7/1, where 7 is the Lorentz factor. 
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To close the hydrodynamical equations we need to 
know the equation of state of the medium. We follow 
Ref.[2lJ where the deconfined phase at high tempera- 
ture is an ideal gas of massless u, d quarks, 150 MeV 
massed s quarks and gluons, and the hadron phase at 
low temperature is an ideal gas of all known hadrons 
and resonences with mass up to 2 GeV[22|. There is a 
first order phase transition between these two phases. 
In the mixed phase, the Maxwell construction is used. 
The mean field repulsion parameter and the bag param- 
eter are chosen as if =450 MeV fm 3 and B x / 4 =236 
MeV to obtain the critical temperature T c — 165 MeV 
at vanishing baryon number density. 

We assume that the medium maintains chemical and 
thermal equilibrium until the energy density of the sys- 
tem drops to a value of 60 MeV/fm 3 , when the hadrons 
decouple and their momentum distributions are fixed. 
We have tested that a slight modification on the freeze- 
out condition will not change our conclusions. 

At SPS energy, we follow Ref.[23j to set the initial 
energy density and baryon number density with the 
wounded nucleon number density, 

n wn (b,x t ) = T A {x t )(l~e-°™ T ^ 3 '-^) 

+ T B (x t - b)(l - e- aNNTA ^) , (3) 

where T A{B) (x t ) = ^^dz p A{B) (x u z) is the thick- 
ness function of nuclear A(B) with nuclear density pro- 
file pa(B)- We choose Woods-Saxon profile, p(r) = 
e (r-H°/; +1 ■ The parameters R and £ for 208 Pb and 
197 Au are from Ref . |24| . At RHIC energy, we follow 
Refs.|25l |26| and assign 75% of the entropy density to 
soft contribution which is proportional to the wounded 
nucleon number density and 25% to hard contribution 
which is proportional to the binary collision number 
density n fcc (6, x t ) = T A (x t )T B (x t ~ b)a NN . The baryon 
density is obtained by adjusting the entropy per baryon 
to 250 26] . The nucleon- nucleon inelastic cross section 
<jnn is 32 mb at SPS and increases to 41 mb at RHIC. 

In peripheral collisions or in the peripheral region of 
central collisions, the system can not reach thermaliza- 
tion even if the energy density is larger than the critical 
value for phase transition. In order to describe this ef- 
fect, we incorporated a cut in the initial condition. We 
assume that the medium will hadronize and decouple 
at the initial time if the local entropy density is less 
than a critical value s c . This critical value will be an 
adjustable parameter of our approach, and will be fixed 
by fitting the J/ip anomalous suppression data at SPS. 
The same value will be used in the calculation of the 
hydrodynamics at RHIC energy. 

We use the well tested RHLLE algorithm^ HU to 
solve the hydrodynamical equations numerically. The 
simple first order operator splitting method is used to 
extrapolate the original one dimensional RHLLE algo- 
rithm to two dimensions. 



III. DISSOCIATION CROSS-SECTIONS 

In the QGP phase, we consider gluon dissociation 
process (g + "J" — ► c + c) only. We will use \& as the 
shorthand notation of charmonium in this Letter. The 



gluon- J/ip dissociation cross section can be obtained 
from the perturbative calculation with non-relativistic 
and Coulomb potential approximation for the cc system 
in the vacuum 291, 



'(uj) =A 



(^-1) 3/2 
(w/e,/,) 5 

-1/2 



(4) 



with A Q = (2 n 7r/27)(TO 3 e^)- 1/2 , where u 
V%Vg^l m ^ — ( s — m^)/(2m*) is the gluon energy in 
the rest frame of \&, = 2m d — Tn^ the binding en- 
ergy of , 3>, m c = m B = 1.87GeV and = 3.1GeV 
are, respectively, the charm quark mass and J/ip mass. 
In a similar way, the gluon dissociation cross section for 
130| and x c @ are 



M = 16A Q 



( W /^'-i) 3/2 (^'-3) 2 



4A (u/e Xc - l) 1 / 2 
(9(^/e x J 2 -20(^/6 Xc ) + 12) 



(5) 



The masses of ip' and \c are 3.7 GeV and 3.5 GeV, re- 
spectively. The J/ip from the feed-down of is about 
10% of the total final J/ip's in pp collisions. For sim- 
plicity, we neglect the contribution of ip' to J/ip in our 
calculation. Therefore, 40% of the final state J/^'s 
are from the feed-down of Xc and others are created 
directly 0- 

We have neglected the medium effect on the dissoci- 
ation cross sections. In relativistic heavy ion collisions, 
the medium produced is quite finite and evolves very 
fast, perhaps the medium effect is not so important as 
that in a static and infinite medium. 

In the hadron phase, the most populated hadrons are 
pions. There are many effective models that can cal- 
culate the inelastic cross sections between charmonium 
and hadrons[3l|. For J/tp, the dissociation cross sec- 
tion is about a few mb which is comparable to the gluon 
dissociation cross section. If the QGP is created in rel- 
ativistic heavy ion collisions, the thermalized hadron 
phase is at the later time of the evolution, and at that 
time the number density of hadrons is very small com- 
pared with the number density of gluons at the early 
time. Although the hadron dissociation may be impor- 
tant to the consideration of \c and ip 1 suppression, we 
neglect the contribution of hadrons in the calculation 
for J/ip. 

IV. CHARMONIUM TRANSPORT 

We define the full phase space distribution function 
of charmonium at a global time t as (x, p, t) . The 
transport equation describing the evolution of is 

P^df, h (x, p, t) = -C% (x, p, t) fm (x, p, t) . (6) 

The l.h.s is the drift term and the r.h.s is the dissocia- 
tion term due to the inelastic collisions between char- 
monium and constituents of the medium. The coeffi- 
cient C* of the dissociation term in Eq. © is 



C B (x,p,t) = i 



d 3 k 



{2nf2E k 



aZ(s)4F g *f g (x,k,t) (7) 
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with Ei. = 



i 2 + k 2 and flux factor F, 



9* 



{p^Pgn) 2 ~ where f g is the Lorentz invariant 

distribution function of the constituents, and c*(s) is 
the dissociation cross section for the charmonium in the 
medium. It is obvious that is Lorentz invariant. 

Considering the boost invariant initial condition at 
the proper time To, it is convenient to use momentum 
rapidity y = 1/2 In {{E +p z )/(E — p z )) and space-time 
rapidity r\ = 1/2 In ((t + z)/(t — z)) as longitudinal 
variables. With the assumption that the charmonium 
is produced at the point t = and z = due to the 
large Lorentz contraction of nucleus and short duration 
time of the initial hard processes in high energy heavy 
ion collisions, the charmonium momentum rapidity is 
equal to its space-time rapidity, y = rj. The distribution 
function can then be rewritten as 

h(xuV,Puy,T) = —fo(x t ,Pt,y,r)5(ri-y) . (8) 

T 

After integration over rj, the l.h.s. of the transport 
equation |JBJ is simplified as 

— ( m t-7^+Pt-J^)fe(xi,Pt,y,T), (9) 
r or oxt 

and the r.h.s becomes 

- — C%(x t ,y,pt,y,T)U(x t ,p t ,y,T) . (10) 

T 

In Bjorken's hydrodynamics, the longitudinal ra- 
pidity of medium flow equals its space-time rapidity. 
Therefore, ^'s rapidity is zero in the longitudinal rest 
frame of the medium flow. Due to the Lorentz invari- 
ance of C*, we can calculate it in this rest frame. Fur- 
thermore, we assume that the temperature and baryon 
chemical potential in the rest frame is rj independent, 
since the rapidity region of the detected 'S's is very 
narrow and near the central rapidity region. With this 
assumption, C* becomes independent of the "J rapid- 
ity, the transport equation can then be simplified as 

( m *^: + Pt- TjL)fa(xt,pt,y,T) 

= -C x l{x u p U T)f^{x t ,pt,y,T) . (11) 
After integration over the rapidity y, it is reduced to 

( TO *J~ + ft' -Jj^)h(%t,Pt,T) 

= -CD(xt,Pt,T)fa{x t ,p t ,T) . (12) 

With Cooper-Frye formula |32^. the number of 'J's at 
proper time r is 



N 9 (t) 



1 



(2tt) 3 

h(xt,v,Pt,v,T) 
l 



d 2 p t d 2 Xt dy drj r m t cosh(y — rj) 



(2tt) 3 



d p t d xtTom t f*(xt,Pt,T) . (13) 



Since the normal suppression of \E' has ceased before 
the starting time of the medium evolution, the number 
of 'S's at tq is 



N*(t ) = / fN(xt,Pt\b)d 2 x t d 2 pt , 



(14) 



where /at is the ^ transverse distribution function at 
impact parameter b after the normal suppression^^, 



f N (x t ,p t \b) 



'NN 



dz A dz B p A (x tl z A )p B (x t - b, z B ) 



e — &abs (TA(xt,ZA, + 0o)+TB(x t ~b, — 0O,ZB)) 

* e -pV{pl) 



(pi 



with averaged transverse momentum square | 



(15) 



(p 2 t ) (b,x t ,z A ,z B ) = (p 2 ) NN + (16) 
a gNPo 1 (r A (x t , -oo, z A ) + T B (x t - b, z B , oo)J , 

where the thickness function T is defined as 
T(xt, zi, z%) — dz p(x~t, z). The constant a a bs is usu- 
ally adjusted to the data from pA collisions where ^'s 
experience only normal suppression. From Ref . |33j| . we 
get the latest value, <r a b s =4.3 mb for J/ip and \ c at 
SPS energy. At RHIC energy, the d-Au experiments 
show that there is litte normal suppression |34|, we take 
Cabs = in our calculation. The p t broadening [35| ef- 
fect induced by gluon-nucleon elastic scattering in the 
initial state is reflected in the second term on the r.h.s 
of Eq. (|16|l . The latest estimate of the constant a g ^ 
is 0.076 GeV 2 /c 2 fm- 1 at SPS energyj^, we assume 
that the value is not changed at RHIC energy and the 
same for J/tp and \c- The averaged transverse mo- 
mentum square (p 2 ) MM in nucleon-nucleon collisions is 
1.15 GeV 2 /c 2 Hg at SPS and about 2.70 GeV 2 /c 2 |13 
at RHIC. Again we use the same value for J/tp and \ c . 

After the comparison of Eq. (|13|) with Eq. Ijl4|l , we 
get the initial condition at To for Eq. I|12fl . 

fn,(x t ,p t ,To) = (2ir) 3 f N (xt,pt\b)/(T m t ) . (17) 

The transverse transport equation l|12f) can be solved 
analytically with the result |l8( 



/*(£ t ,pl,T) = e 



- LJ <1t'C% (x t -v t (r-r'),p*t ,r') / m t 



X h(x t ~ Vt(T - T Q ),pl,T ) 



(18) 



It should be noted that the leakage term, that is the 
second term on the l.h.s. of Eq. l|12fl . has been shown 
in Ref. 01 to be important in the consideration of char- 
monium p t distribution. 

When we know the time evolution of the produced 
medium and the ^ dissociation cross section in the 
medium, we can calculate the number of and its p t 
distribution at freeze-out time t/ with Eq. (|13|) . 



V. NUMERICAL RESULTS 

We first fix the parameter s c in the initialization 
of hydrodynamics by fitting the J/ip suppression data 
from SPS. The best fit is achieved with s c — 31.7 fm -3 . 
Figd (A) is our calculation compared with the NA50 
data[33|. The thin dashed and solid lines indicate, re- 
spectively, the theoretical calculation of normal sup- 
pression with Drell-Yan cross section rescaled to MRS 
43 and GRV LO parton distribution function |33l|. The 
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FIG. 1: The J/tp suppression and (p? ) as functions of cen- 
trality at SPS energy. 



FIG. 2: The J/tp suppression and (p? ) as functions of cen- 
trality at RHIC energy. 



thick dashed and solid lines are the corresponding re- 
sults with anomalous suppression. The transverse en- 
ergy E t = 0.27AN p {b) GeV 37] is the measure of cen- 
trality with N p (b) being the average number of partici- 
pants in the collisions with impact parameter b. In very 
peripheral collisions with Et < 23 GeV, the maximal 
initial entropy density is less than the critical value s c , 
the medium will decouple at the initial time, and there 
is only normal suppression. It is obvious that due to the 
missing of E t fluctuation in our calculation, the anoma- 
lous suppression starts suddenly at E t — 23 GeV. With 
increasing centrality, the volume and the life-time of the 
produced QGP will increase, and the anomalous sup- 
pression becomes more and more important. A good 
fit with the suppression data has been achieved in the 
whole E t range except in most central collisions where 
the fluctuations was shown to be important [37|. 

After fitting the parameter s c , we can calculate ob- 
servables related to the pt broadening effect. In FigQ] 
(B), we show the (pf) as a function of E t . The 
dashed and solid lines are, respectively, J/tp (p^) cal- 
culated without and with anomalous suppression. The 
later agrees well with the experimental results from 
NA50 36]. With increasing centrality, (pf) increases 
steadily till E t — 90 GeV. After that, (p^) becomes 
saturated and finally drops slightly in most central col- 
lisions. There are two competing mechanisms that af- 
fect J/tp (p 2 t ): The J/tp with large Pt is mostly pro- 
duced in central collisions according to the Cronin ef- 
fect. Since the matter produced in central collisions is 
denser and hotter than that in peripheral collisions, the 
J/tp with large pt has more chance to be absorbed by 



the QGP; On the other hand, the J/tp with large pt 
has more chance to escape the anomalous suppression 
region, when its duration time in the QGP is shorter 
than the anomalous suppression time. This is the leak- 
age effect initially indicated by Matsui and S atzp and 
studied in detail in the transport approach in |l8( . The 
former mechanism suppresses the (p? ), while the later 
enhances the {pf)- 

With the fixed parameter s c , we calculate now the 
J/tp suppression and {pf ) at RHIC energy. FigEl (A) 
shows that our result on J/tp suppression agrees with 
the low statistical data from PHENIX[3g. The sup- 
pression at RHIC is much stronger than that at SPS, 
but does not approach zero. Since there is still no trans- 
verse momentum data from RHIC, we displayed only 
the theoretical results in Fig|2] (B). The behavior of 
the (pi ) at RHIC is similar to the corresponding result 
at SPS energy. However, it is obvious that the leak- 
age effect prevails even in very central collisions, and 
there is no clear decrease of (p?), very different from 
the predictions 39, 40] calculated with p t independent 
dissociations and without considering the leakage effect. 

Finally we consider the elliptical flow defined by 



'-'2 




(19) 



In high-energy nuclear collisions, the initial spatial 
anisotropy is transferred to final momentum anisotropy 
via partonic and/or hadronic interactions. Interesting 
information about the early stage of the hot and dense 
medium may be obtained from the study of this az- 
imuthal anisotropy 0, E2> HE H 03 • The measured 
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FIG. 3: The elliptical flow V2 of J/tp as a function of cen- 
trality at SPS and RHIC. 



values of anisotropy v 2 of hadrons with light quarks 
(u,d,s) at RHIC are positive and large compared with 
lower energy results, and show a clear hadron mass 
dependence in the low transverse momentum region, 
pt < 2 — 3 GeV/c 0,^3- Hydrodynamic calculations 
have been used to explain the observed collective expan- 
sion behavior 45]. The intrinsic charm mass is large 
compared with the initial temperature that might be 
reached at RHIC |4q. Therefore, any collective motion 
of charmed hadrons will be a useful tool for studying: 
(i) charm production mechanism in nuclear collisions: 
via direct earl y pQCD -typ e creation ^ij or via later 
coalescence [5(j, 13, EH l53l| : (ii) early thermalization of 
light flavors (u,d,s) in the high-energy nuclear collisions. 
In the case of thermalization, if coalescence is the dom- 
inant process for charmed hadron production, v 2 could 
be largejSl], while it is predicated that the direct pQCD 
produced charmed hadrons carry small or zero v 2 . How- 
ever, the leakage effect, discussed above, will lead to fi- 
nite v 2 for J/tp in non-central collisions [a, [EH- This can 
be understood transparently: The anormalous suppres- 
sion of the J I tp depends on the length that the particle 
travels throught the medium. Hence more suppression 
is expected in the out-of-plane direction (y-dircction) 
which leads to finte value of v 2 with positive sign. 

In FigElthe solid and dashed lines are our results of 
J/tp v 2 at fixed transverse momentum p t — 3 Gev/c at 
SPS and RHIC. It is clear that v 2 is not zero but finite in 
mid-central collisions. For very central and very periph- 
eral collisions v 2 approaches to zero due to symmetric 
anomalous absorption and no anomalous absorption, 
respectively. The values of v 2 at RHIC are found to be 
larger than that at SPS due to the fact that the colored 
medium created at RHIC has larger volume and longer 
life-time which lead to the stronger leakage effect. If the 
J ftp's are created at hadronization, their minimum-bias 
v 2 has been evaluated in coalescence model with the 
assumption of complete thermalization of charm quark 
with the medium and the assumption of the same max- 
imum v 2 (about 9%) for charm and light quarks 51j. 
The result is shown in Fig. 21 as a function of pt and 
compared with our calculation at fixed impact param- 
eter b = 7.8 fm, at which the v 2 is found to have the 
maximum. The minimum-bias v 2 which is the average 
of v 2 over the impact parameter should be less than the 
maximum. It is clear to see that the maximal v 2 calcu- 
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FIG. 4: The elliptical flow of J/tp as a function of pt at 
RHIC energy. The solid line is the maximal V2 with impact 
parameter b=7.8 fm calculated in the frame of J/tp trans- 
port, and the dashed line is the minimum-bias V2 (scaled by 
a factor of 0.1) of the coalescence model with the assump- 
tion of complete charm quark thermalization. 



lated in the frame of J/tp transport is less than 10% of 
that calculated in the approach with full charm quark 
thermalization 15 J. 



VI. CONCLUSIONS 

By combining the hydrodynamic equations for the 
QGP evolution and the transport equation for the pri- 
mordially produced J/tp in the QGP, and considering 
the anomalous charmonium suppression induced by the 
gluon dissociation process, we calculated the J/tp sup- 
pression, averaged transverse momentum square, and 
elliptic flow at SPS and RHIC energies. It is found 
that the leakage effect reflected in the transport equa- 
tion and the p t dependence of the charmonium dissocia- 
tion cross sections play an important rule in explaining 
the experimental p t broadening effect at SPS. 

In our transport approach, the main driving force of 
the elliptic flow is the leakage effect, the J ftp's mov- 
ing in the direction of the long axis of the anisotropic 
medium are strongly suppressed, but those moving in 
the direction of the short axis are easy to escape the 
medium. The inelastic interaction with the medium 
will transfer collective flow to the charmonium, but the 
effect is small. Since the primordially produced Jf tp is 
not thermalized in our transport approach, the calcu- 
lated v 2 is the low limit of the J/tp elliptic flow. Any 
calculation with partial or full charm quark or char- 
monium thermalization should be larger than this low 
limit. From the comparison of theoretical calculations 
in different models with the experimental data, the el- 
liptic flow may be a useful tool to determine which 
mechanism is the dominant one of J/tp production, ini- 
tial pQCD production, coalescence at hadronization, or 
the mixture of the both. 

The possible gain term on the right hand side of the 
transport equation due to the recombination of c and c 
in the parton stage will probably enhance the charmo- 
nium production at RHIC energy and affect (j% ) and 
elliptic flow. These effects are now under consideration. 
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